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X. INTRODUCTION 


This report describes technical findings of work supported 
by contract no. 956525 for the period, July 1, 1984 to Sept. 31, 
1984. Details of the findings appear in the Appendix. In this 
section, we briefly summarize these findings and indicate their 
practical implications for salat — cell design. 

A controversy exists about several of the key parameters 
describing the heavily doped regions of Si solar cells. This 
controversy relates particularly to experimentally determined 
energy-gap narrowing and to the minority-carrier diffusivity and 
mobility. To explain the origin of this discrepancy, 1st us 
denote the energy-gap <or bandgap) narrowing by aE and the 
minor! ty-carri er diffusivity and diffusivity by D and p. We omit 
subscripts for brevity. 

1.1 Some Equations for Heavily Doped Silicon 

All experimental determinations of AE, except those 
employing photol umi nesence or optical absorption, depend on 
measurement of mi nori ty-carri er current Jin the heavily doped 
region. If we assume n-type heavily doped Si for illustration, 
this current is given by 

J => PPVE V - eDvP < 1 ) 



is the quasi -f J eld acting on the minority holes. Note that the 
field acting on the holes differs from that acting on the 
electrons and differs also from the electric field <the gradient 
of the electrostatic potential). Eq. (2) follows from the 
definition of the energy gap and from the condition, from the 
quantum theory of solids, that a particle at a band edge has only 
potential energy (no kinetic energy) . 

For low-injection conditions, which hold always in a heavily 
doped region of a Si solar cell, 

P(r) =Cn 2 /N (r ) 3 expCF - F 3 (3) 

i e n p 

where the argument of the exponential in Eq. (3) is the 
difference between the electron and hole electrochemical 
potentials (or quasi-Fermi levels) normalised by kT/e (or kT) . 

The effective intrinsic density is given by 


n 2 5 NR » n 2 expCAE'3 F. .-ME' - £')3 + expCaE'3 (4) 

le o o x 1/2 f c 

where the primes mean normalisation by kT and where F^ 0 denotes 
the Fermi-Dirac integral of order 1/2. 

A simple, nearly tutorial, development of the facts stated 
above appears in Ref. 1 for this Section. 

1.2 Measured Recombination Current 


Now we are in a position to examine the measured 


recombination current J of a heavily doped region; we continue 
assuming n-type Si far illustration. By integrating the hole 
continuity equation across the quasineutral portion of the 
heavily eloped region, we obtain 

J/b - ,f(AP/T)dx + J (surface) (51 

in which AP is the hole density in excess of the equilibrium 
density, t is the position-dependent hole lifetime, and 
J (surface) is the hole recombination current at the surface. Eq. 
(5) displays the hole recombination current of the n+ region as a 
sum of volume and surface recombination. 

One can measure J by EBIC, by dark current in presence of 
applied forward voltage, or by selectively illuminating the 
quasineutral n+ region and measuring the short-circui t current. 
However the measurement is made, the problem in using J to 
determine AE lies in separating the volume and surface 
recombination components of (5) and connecting then to AE via 
Eqs. (1) to (5). 

1.3 A Key Assumption and a Questionable Tradition 

From Eqs. (1), (3) and (4), we see that D, p, and AE 

interplay to determine J. Thus, in principle, one must determine 
D and p (linked by D/p ® kT/e) separately from AE. 

It has been the tradition in all but recent work by Arnost 
Nuegroschel and me on heavily doped Si to avoid determination of 


D and ,u. Rather, it has been assumed that the majority carrier 
and minority carrier D and M are equal. This assumption 
underlies all models -for AE that are widely used, beginning with 
the first experimental determination by Slotboom and deQraaf . In 
particular, this includes the models now in computer programs at 
3 PL. 

To my knowledge, there is no reason, experimental or 
theoretical, for making this assumption. It is a convenient 
assumption because the majority carrier mobility is easily 
determined experimentally. Minority-carrier mobility, prior to 
our recent work, was only measured for dopant concentrations up 
to lO^/cm 3 . From a theoretical viewpoint, Bennett (1983) of the 
National Bureau of Standards has proposed that the majority 
carrier mobility is less than the minority carrier mobility by 
reason of the sum rule of quantum theory. Kane (1984) of Bell 
Laoorator i es , in a preprint sent to Prof. Neugroschel , has 
indicated his belief that the majority carrier mobility exceeds 
the minority carrier mobility. His reasoning relies on the simple 
observation that the kinetic egergy of the majority carriers able 
to conduct greatly exceeds the kinetic energy of the minority 
carriers. Hence, in a quasi-cl assical picture, the majority 
carriers spend less time in the vicinity of the force field of a 
scattering center such as an impurity ion. Consequently they are 
scattered less than are minority carriers and their mobility (JJ 
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er scattering' /m * ) ma V be B >*l? ec tsd to be larger. 

1.4 Suggestions -from our Recent Experiments 

In our work on the experimental determination of AE and p 
and D (of the minority carriers) , we avoided the assumption of 
equal majority and minority carriers by exploring the temperature 
dependence of J in a small range of temperature. This led to Ref. 
2 (below) on AE and to Ref. 3 (below) on and D. In these 
methods dClog Jl/dT^is explored. The activation energy thus 
revealed yields aE (1982); use of this AE together with the value 
of J then yields the minori ty-carri er D and p (1983). This work 
in 1983 suggested evidence for minori ty~carrier mobility that was 
about one order of magnitude below the majority carrier mobility . 
We suggested a model consistent with this evidence, in which 
minority carriers made transitions to trap levels energetically 
near the edge of the minority carrier band edge (mobility edge) 
and subsequently were released. This removal from the conducting 
state, for a time characterised in the Appendix by T ^ ra ppi n gJ 
decreases the minority carrier mobility. The model of Kane also 
leads to this same result. Either model may explain the 
experimental evidence. 

1.5 Contents of the Appendix 

This model is explored further in the Appendix. There we 
excite by short wavelength illumination and demonstrate 
activation behavior of the mi nori ty-carri er mobility and 



diff usi vity, Short-circui t current is the measured variable. In 
our interpretation , the activation energy is determined -from the 
measurementsto be of the order of h.T , which is expected for 
bandtail states uut does not rule out an acceptor level as the 
dominant trap. 

We argue further that the short circuit current times 
temperature to the mth power, that is 

I sc T m vs i/T, 

shows activation behavior. Thus, based on our experiments, the 
assuintion that /J (majority) ~ M (minority) is always wrong. This 
is implied in Eq. <B) of the Appendix. 

1.6 Conclusion 

From a technical standpoint, there is no conclusion or 
anything approaching a conclusion at present. We have suggested 
evidence for our findings for AE which are inextricably linked to 
our findings for the rninori ty-carrier mobility and diff usi vity. 
The Appendix of this report extends that evidence. But we 
recognize that concentration on activation energies, for the 
problem under study, relies on values of slopes deriving from 
current measured as a function of 1/T. A small error in the 
slope can produce a large error in the measured aE and D and ,4. 

We continue work toward decreasing this error. 

Other workers have dogmatically insisted that values of AE 
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of Slotboom and deGraaf and subsequent workers, which are 
markedly lower that our suggested values, are in -fact correct. 
This is done even though there is no reason to believe the 
equality of majority and minority-carrier mobilities for 
concentrati ons well above iO^/em 8 *. This equality is equivalent 
to believing in the lower values of energy-gap narrowing. 

1,7 Practical Consequences 

The performance parameter of silicon solar cells that 
probably links most closely with aE is the open-circuit voltage. 
For existing solar cells, use of our values of AE (vs impurity 
concentration) in any of the common computer codes will probably 
yield the measured open-circuit voltage, to a good approximation, 
provided one uses our values of /-t and D. On the other hand, use 
of the more widely accepted values of AE together with the 
assumptl on that Mmajority) - ji(minority) will probably also 
yield the measured open-circuit voltage. The reason for this one 
can see by inspection of the equations given above. Open-circuit 
voltage is determined by J, primarily? and J derives from P, 
which is determined by AE, and from D and !i, together with the 
quasi -field of Eq. (1) which depends on AE (Eq. 2). 

Thus this simple test will not discriminate between the 
adequacy of the models. Better tests may include: 

(a) sensitivity of the open-circuit voltage to the surface 
recombination velocity S. 


(b) temperature dependence of the open-circuit voltage. 

It is not clear whether the spectral response will provide an 
adequate test for the models because of the connections among aE, 
P and D that underlie the spectral response. The temperature 
dependence of the spectral response may be more revealing. 

Indeed, the Appendix explores the temperature dependence of the 
short-wavelength current; the tentative result there favors the 
data of Ref. 3 rather than the assumption of the equality of the 
majority- and minori ty-carri er mobilities. 

Because many solar cell designs of the present depend for 
their quality on the surface recombination velocity, much more 
work is needed aiming to answer the questions about AE, p and D 
raised here. It is probable that the values for AE, P and D 
determined in Refs. 2 and 3 are inaccuratP, for reasons 
discussed. So also are the values of these parameters widely in 
use (including in JPL computer codes). Until more accurate 
values become available, computer-assisted evaluation of new 
designs will be misleading, as will the determination of values 
of other parameters such as the surface recombination velocity 
bordering a highly doped emitter region. 
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APPENDIX 


TRAP CONTROLLED MINORITY-CARRIER MOBILITY IN HEAVILY DOPED SILICON 

A, Neugroschel, F.A. Llndholm 
Department of Electrical Engineering, 

University of Florida, Gainesville, FL 32611 
and C.T. Sah 

Department of Electrical Engineering, 

University of Illinois, Urbana, IL 61808 


Recently [1], we presented experimental evidence for low diffuslvity and 
mobility of minority holes in highly arsenic doped (concentrations 
~ 10 20 cm" 3 ) n-type Si, This evidence suggests that the minority-carrier 
dlffusivity (and mobility) may be about an order of magnitude smaller than the 
majority-carrier hole dlffusivity in comparably doped p-type Si. As an 
explanation, we suggested a simple transport model that emphasized trapping by 
localized tail states in the minority band. The commonly used assumption of 
equal minority and majority-carrier 1 mobilities and the neglect of the tail- 
states effects on the minority-carrier mobility was also questioned in a 
review article by Abrams et al . [2], Here we report qualitative evidence that 
minority-carrier diffusion in n* Si may be trap limited. 

Figure 1 illustrates the band structure, including the tail states, for 
n + Si. As was already discussed in [1], for low-level injection, the minority 
holes occupy energy levels near the mobility edge Ey. In the model proposed, 
the holes from the extended states in the valence band can be captured by the 
localized tail-states for some mean time t^pp^g, and then released back 
into the valence band. This process, shown by arrows in the Fig. 1, will 
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decrease the hole mobility. If the scattering rate of holes Inside the band Is 
comparable with the hole trapping rate, It Is useful to generalize the 
Boltzmann equation to include band-bound transitions. This generalization, 
discussed, for example, by Smith, Janak, and Adler [3], and worked out in 
detail by Sah and Lindholm [4], involves approximating the collision integral, 
for small departures from equilibrium, by the sum of two terms of the form, 

i 

(f - f 0 )/T. Here the first term Is the customary relaxation time 
approximation for the nonequilibrium distribution function in which r scatt Is 
the scattering time that characterizes intraband transitions. In the second 
term, Tt rap p-| n g is mean trapping time for localized states near the valence- 
band edge. Thus 

"p ■ e W l, p > 1/T coll " >/’,c.tt + ‘/’trapping • 0) 

We inve'iti;!f.ce the temperature dependence of u p by estimating T trapp1ng C4]: 

'/‘trapping a e p “ atI " ex P [ -( e t * E V )/kT] < 2 > 

where A Is a temperature Independent constant. In (2) we have assumed, for 
simplicity, trapping at one shallow level (close to E v ) in the tail band or at 
an acceptor level in the n + region. Such an acceptor level can come from the 
dopant acceptor (boron in our case) In the p-type substrate of a diffused or 

ion implanted n + /p diode. For a parabolic valence band and no degeneracy, 

a 

m = 2 [4]. However, for A band-tail distorted band, the value of m is not known 
and must be determined experimentally, as is attempted below. At low and high 
temperatures, respectively, (1) reduces to 


- 2 - 



1/,T scatt ’ l/,T co1 1 1 hi gh T a ^trapping* 


(3) 


1/T coll|low T a 
If 1/T trapping Is dominant, (l)-(3) suggest 

trapping) “ C(e/n.*)/ { AT m ) ] axp(E A /kr) 


(4) 


where E/\ = Ey - Ey > 0 is the activation energy of the bound-state level Ey. 

The physical interpretation of (1) - (4) is as follows. In heavily doped 
n-type Si, holes reside mainly In the delocalized states of the valence 

band. Only a small fraction of the holes is trapped in the localized states 
of the narrow band tail. At low temperatures, the holes in the tail states 
are frozen there. Thus u p is determined entirely by the intraband 

scattering: l/x coll ~ 1/ T scatf As temperature Increases, trapping 

transitions become important (see Eq. 2) and u p decreases (see Eq. 1), which 
allows E^ to be deter; ained using (4). 

Experiments were done to provide qualitative support for the hole- 

trapping model. The experiments involved a measurement of the temperature 
dependence of the photocurrent response of both n + /p and p + /n photodiodes. 

The wavelengths of the incident light used create electron-hole pairs almost 
entirely in the heavily doped n + or p + region. Transport of photogenerated 
minority carriers to the collecting p/n junction and the consequent jhort- 

circuit current I sc depends on the minority-carrier mobility jj and diffusivity 
D [5]. Hence I sc and its temperature dependence supplied a vehicle for 
studying u and D. 

The devices for which we report findings here were fabricated by arsenic 
implantation into 5 n-cm p-type Si substrates followed by 1200°C anneal for 

30 minutes. The resulting n + layer was 1.2 jim deep and the sheet resistance 
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was 10 n/square. The As concentration was ~ 10 30 cm“ 3 and nearly Independent 
of position x over about 0.6 pm below the surface (x *» 0); for x > 0.6 pm it 
slowly decreases becoming about 3 x lO^cm -3 at 1 pm. The surface was 
unpassivated, covered only by a thin {~ 10 A) native oxide layer. Other 
devices with both unpassivated and S10 2 passivated surfaces vere also 
studied. The results are similar to those to be presented here. The metal 

I 

contacts covered less than 10X of the area of the front surface. 

The Isq{T) dependence was measured in the wavelength range X = 0.38 pm to 
0.4 pm for which the electron-hole generation rate follows G(x) = G(0)exp{-ax) 
[5], where a > 10 5 cm*** [6], Hence contributions of G(x) to Ij C originating 
In the p/n junction space-charge region and in the p-type substrate are 
negligible. Thus [5] for (a L ) 2 » 1 

r 

AqF ( 1 - R) aD p + S p 

! SC ** air TT~ . 

p D ({“) + C0S ^ 

p p p 

Here A is the device area, F is the illumination flux density, R is the 
reflection coefficient, Sp is the effective hole surface recombination 
velocity, W - 1.3 pm is the thickness of the quasineutral n + layer, 
Lp = (OpTp) 1 ^ 2 is the hole diffusion length, x p is the hole lifetime, and D p 
is the hole diffusivity. 

The analysis of (5) is complicated if one wishes to obtain the magnitude 
of the diffusivity D p or mobility p p = 0 p (kT/q) from the measured I sc . In 
this note we are interested, however, only in qualitative trends and try to 
present a simple picture consistent with the data. For this purpose we 
consider now the temperature dependence of I^. The temperature dependences 
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of a nd R [7] and tp * t a {Auger lifetime) [8] are very small. This leaves 
Sp and Dp as the only temperature dependent parameters in (5). We neglect the 
temperature dependence of Sp, for a moment, later justifying this assumption 
based on experimental results and on theoretical predictions. 

Recent measurements indicate that Lp < 1 pm in the n + silicon doped at 
about 10^0 cm“ 2 [9]. Since W/L “ 1 (where W « 0.5 pm is the width of the 

r 

highly doped portion of the emitter at the surface) we can approximate (5) as: 
ctD + S 

*SC “ B aLS p +(Dp/L p ) j {6) 

where B = AqF(l-R)/LpSinh(W/Lp) is almost temperature independent. Consider 
(6) for three special cases. First, if Sp >> aDp and Sp » D p /Lp , l sc « B/a. 
In this case I S q f- f(T). This contradicts our experimental result that I S q 
increases with T (Fig. 2), Second, if- S p « D p /L p and S p « aD p , I sc « BL p = 
B(TpDp) 1 / 2 , i.e. I sc * D p i//2 - Third, if (Dp/L p ) < S p < aD p (for our device 
2x10'^ cm/sec < S p < 10 5 cm/sec, with D p ~ 0.2 cm 2 /sec [1], a > 10 5 cm -1 , 

Lp = 10" 4 cm)* then I sc = B(D p /S p ) . We assume now that the third case 

prevails in our devices, although the analysis below will also accomodate the 
second case and thus the range 0 < S p < 10 5 cm/s. Lastly, the assumption 

(aLp) 2 » 1 used to derive (5) is also valid, since (aL) ~ 10. Hence, the 

analysis is self-consistent. 

To obtain simple picture, we try the assumption, m=0, and find from 
(4) and (6) that 


I sc /T = C exp(E A /kT) 


( 7 ) 
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moves Epp closer to Ey (Fig. 1) which increases the occupation of the tail 
states closer to Ey and E^« (Ey - Ey)/kT is expected to decrease, as 
observed. In the temperature range from about 70 K to about 160 K, where the 
mobility activation behavior is observed, the variation of the energy gap with 
T is very small [11], which justifies our assumption of nearly temperature 
independent ct and R [7]. 

The effect of the temperature variation of the effective surface 
recombination velocity Sp(T) on I SC (T), which was neglected, is now made 
credible by the following considerations. First, if Sp ( T ) is important, E^ is 
expected to be about Eq/ 2 « 0.55 eV in a Shockley-Read-Hall model even though 

the surface has distributed states in the energy gap [12]. Second, Sp at the 

heavily-doped n + surface is likely dominated by temperature insensitive Auger 
effects [12]. Third, results similar to those shown in Fig. 2 were also 
obtained from the devices for which passivation by Si0 2 sharply decreased Sp. 

The fact that I^q is not negligible below ~ 100 K indicates that most of 
the minority holes are in the extended band states at T ~ 100 K with only a 
small fraction trapped at the tail states and immobilized. For the 
excitations used to measure Ig^T) and also for low-level hole injection in 
the dark u c ed in [1], the hole concentration in the n + region is very small. 
We can then conjecture that the penetration of the localized tail states into 

the energy gap is only a very small fraction of the energy gap. This 

supposition, based on our data, agrees with recent theoretical calculations of 
the band structure of heavily doped Si [13], 

It follows directly from our trapping model that ii(minority) < 
u(majority) could apply also in the absence of the tail states because shallow 
impurity levels can act as traps near the minority-band edge (Fig. 1). This 
suggests that u(minority) may strongly depend on compensation. A lack of 


- 7 - 



compensation in their epitaxial p/n diodes may be responsible for the 
experimental observation by Dziewior and Silber [14] that li(minority) - 
p(m.ajority) for concentration in Si < 10* 9 cm" 3 . 

Although the assumptions needed to derive (7) from (5) introduce some 
inaccuracy in the model for the activation energy of minority-carrier mobility 
(diffusivity) , the general conclusions derived by comparing (7) with the. 
experimental trends exhibited in Fig. 2 are anticipated to remain valid. The 
exact value of the activation energy in a single-trap model remains in 
question. When more detailed knowledge concerning the values and the 
temperature dependencies of parameters in (5) become available, one can then 
use data presented here to explore aspects of the minority-carrier band tail. 

The purpose here is less ambitious. We have demonstrated activation 
behavior of the minority-carrier mobility and diffusivity, exhibiting thereby 
an activation energy of the order of kT (for 300 K) as one expects for the 
prominently active tail states. Moreover, our data in Fig. 2, as interpreted 
here, supports the inadequacy of the commonly used assumption of equal 
majority- and minority-carrier mobilities; in turn, this supports 

interpretations, such as that in [1], of larger values of energy-gap narrowing 
than is common in the literature. 

Lastly, we point out that the assumed value, m = 0, and values m < 0 in 
(4) are self-consistent with the experimental results of Fig. 2, but values 
m > 0 are not. However, regardless of the value of m, the I sc T m vs. 1/T plot 
shows an activation behavior, thus the assumption of y(majority) = y(minority) 
is always wrong (see Eq. 8). 

This work was supported in part by National Science Foundation Grant 
ECS-8203091 (A.N*) and the Jet Propulsion Laboratory under Contract 

No. 956525. We acknowledge useful discussions with Peter T. Landsberg. 
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Fig. 1. Qualitative Illustration of the band edges of heavily doped 
n + -silicon. The broken lines show the unperturbed parabolic bands. The 
positions of both the electron and hole quasi-Fermi levels is also 
indicated. The arrows near E v indicate hole capture and emission by the tail ( 
states and by the acceptor level from the p-type substrate. The penetration 
of the tail states into the forbidden gap is assumed to be very small in 
comparison with the bandgap Eq = E c - E v . 

Fig. 2. Normalized short-circuit photocurrent versus 1/T for X = 0.4 pm. The 
illumination density for curve 1 is five times larger compared to that of 
curve 2. 
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956614 

Univ. of Washington 

Annual 9-1-83 
to 8-31-84 

9950-1001 

W08746-83-1 

Solar Energy Res. Inst. 

Qtrly. 7/15/83 
-10/15/83 

9950-909 

956797/02 

SPIRE 

Qtrly. No, 2 

9950-969 

956477 

Univ. of Wise. - Madison 

Final 

9950-1011 

956786 

Westinghouse R&D Center 

3rd Qtrly. 

9950-1010 

956841 

MIT 

3-21 to 5/20/84 9950-1021 

956541 

HR Textron, Inc., Systems 
Eng. Div. 

Extension Final 

9950-1023 

956797 

SPIRE 

Qtrly. No. 3 

9950-1025 

956042 

TRW 

Final Test 

9950-1027 

956384 

Dynamics Research Corp. 

Technical 

9950-1030 

955637 

Dynamics Research Corp. 

Final 

9950-1029 


,-',OGt - 


FV 


Contract No. 

Corporate Source 

Report 

No. 

File No. 

*955637 

Ford Aerospace & Comm. Corp Test 


9950-995 

956722 

Ford Corp. 

FR 


9950-1048 

956885 

RDA Logicon 

FR 


9950-1047 

956503 

System Planning Corp. 

FR 


9950-1051 

955678 

University of Arizona 

FR 


9950-1052 

*956457 

Aerojet Tech Systems Co. 

Final 


9950-974 

956428 

Arizona State University 

FR 


9950-1026 

956525 

Univ. of Florida 

Qtrly. 


9950-1033 

956831-3 

Arco Solar, Inc. 

Qtrly. 


9950-1035 

956064 

SRI International 

Final 


9950-1062 

956042 

TRW 

Final 


9950-1070 

95603B 

Hughes Aircraft Co. 

Final 


9950-1069 

956909 

Earth Satellite Corp. 

Final 


9950-1078 

954349 

Thermo-Electron Corp. 

Final Technical 

9950-1080 

956312 

Mobil Solar Energy Corp. 

Qtrly. 10/1 - 
12/ 31/84 

9950-1036 
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